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Abstract 

We show that several results in seemingly independent areas of num- 
ber theory such as divergent series, summation of arithmetic functions, 
uniform distribution modulo one and summation over prime numbers can 
be unified under a single equation. We introduce the concept of asymp- 
totic equidistribution and apply our method to derive several new results 
in the above areas of topics, especially in the theory of prime numbers. 

1 Introduction 

In 1873 James Clerk Maxwell combined electricity and magnetism using the 
Maxwell's Equations into a single force called electromagnetism and since then, 
physicists have been trying to discover the elusive unified theory that can ex- 
plain all the forces in nature using a common set of equations. Although the 
title of our paper may appear too broad at this stage, we present the analogous 
idea of unifying seemingly unrelated topics in number theory such as the di- 
vergent series, summation of arithmetic functions, uniform distribution modulo 
one and summation involving prime numbers, under one common equation. 

We shall see that several results in the above areas of mathematics which are 
currently treated as independent results can be unified under a single equation 
using summations of the type Y]"i d r f(S r ) where S r = Yl^—i d t . In section 3, 
we present our unified Master Theorem. In section 4, we present the Equivalence 
Theorem which shows that the ratios of the partial sum ^ of a divergent series 
has properties analogous to that of a sequence a n uniformly distributed modulo 
one. In section 5, we introduce the concept of asymptotic equidistribution mod- 
ulo one and prove the sequence of prime numbers is asymptotic equidistributed 
modulo one. In section 6 and 7, we apply the Master Theorem and its varia- 
tions to derive several new results on summations involving prime numbers and 
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derive beautiful asymptotic formula such as: 



lmi —a— > — - + — tan 1 — e Hr+ *i + - + 'r = — - — \G - -^- L ) 

n-s-oo \xf n ' 7 r \ r p r J \ 7 n / 4 V 47T / 

r— 1 

which connects Apery's constant C(3), Catalan constant G, Euler-Mascheroni 
constant 7, exponential constant e, harmonic numner H n , Meissel-Merten con- 
stant M, n th prime number p n , imaginary part of the n th non-trivial zero of 
the Riemann zeta function j n and the fundamental constants of a circle 7r. 
In section 9, we present the generalized the definitions of the Riemann zeta 
function and the Euler constants. All the identities that we have discovered 
in this paper through rigorous methods have also been numerically verified, 
thereby confirming our theory of unification. We have calculated the constants 
that appear in different equations by direct computation using prime numbers 
up to 2 44 . In section 10, we give the details of the verifications of identities as 
well as computations of constants. To conclude our paper, in section 11, we give 
an outline of some of the work that can be done in future in this line of research. 

And with due apologies to Abel, we shall see in this paper divergent series 
are not a useless invention of the devil after all. Divergent series act as the 
bond that unifies different areas of mathematics. 



2 Preliminary lemmas 

In the rest of this paper, a n shall denote a sequence uniformly distributed mod- 
ulo one. For the sake of simplicity, we assume that < a n < 1. If a n > 1 then 
its fractional part {a n } is to be used. 

Lemma 2.1. If a n is uniformly distributed modulo one then for every Riemann 
integrable function f in [0, 1], 



1 - f 1 
lim - YV(a r ) = / 



f(x)dx. (1) 



Proof. Well known. (See g], Page 3) □ 
Lemma 2.2. If f is monotonic and continuous and defined in [l,n] then, 



Y J f{r)= / f(x)dx + 0(\f(n)\ + |/(1)|). 



(2) 



Further if f is continuously decreasing to zero then, 

>(/(»))■ (3) 



£/(r)= / f(x)dx + C f + 0( 

r=l •' 1 



where Cf is a constant that depends only on f. 

Proof. Well known. (See [3], 1.62-1.67, Page 19-20) □ 
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Lemma 2.3. LetS(n) =J2r=i^( r )- V9 is defined in [S (1) , S (n)} andb(x)g(S(x)) 
is monotonic and continuous in [l,n] then, 



^2b(r)g(S(r)) = / g(x)dx + 0(6(n)) 

r=l J S{1) 



r=l 

where 0(S(n)) = \b(n)g(S(n)\ + \b(l)g(S(l))\. 
Using [5] we can express [6] as 

rS(n) 



(4) 



where S(n) = \b(n)g(S(n)\ + \b(l)g(S(l))\ + b'(y)g(S(y))dy. 
Proof. From[2j we have 

S(y) = X>(r) = [ V b(t)dt + 0(\b(y)\ + \b(l)\). 

r=l Jl 

Differentiating S(y) we get S'(y) = b(y) + 0(b'(y)). Let S(y) = x; then, 

S'(y)d y = b(y)dy + 0(b'(y)dy) = dx 

b(y)dy = dx + 0(b'(y)dy) 
g(S(y))b(y)dy = g(x)dx + 0(b' (y)g(S(y))dy) 

b(y)g{S{y))dy = g(x)dx + OU b' (y)g{S{y))dy) (5) 

Again, from[2j we have 

n pn 

^2b(r)g(S(r)) = b(y)g(S(y))dy + 0(6 (n)) (6) 



Y,b{r)g{S{r)) = g(x)dx + 0(S(n)+ b' {y)g(S(y))dy) . 

r=l J\ 

The proof is completed by redefining the error term 6(n) as 

S(n) = \b(n)g(S(n)\ + \b(l)g(S(l))\ + f b' (y)g(S(y))dy. □ 



Notice that [4] is a generalization of [2] If b(n) = 1 then, in the error term 
S(n), the integral vanishes because b'(t) — and so[4]reduces to[2j In this paper, 
we are interested in unifying different areas of number theory rather than on 
developing a general summation formula. So improving the error term S(n) is 
not in scope of our current paper. 

Lemma 2.4. // a n is uniformly distributed modulo m and f(n) is divergent 
then as n — > oo, 

n n 

£a r /(r) = (l + (l))™£/(r) (7) 

r— 1 r— 1 
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Proof. Since a n is uniformly distributed modulo m, (m — a n ) is also uniformly 
distributed modulo m. Hence as n — >• oo, 

n n n n 

E°v/( r ) ~^(m-o r )/(r) = m E /(r) - E a r /(r), 

r— 1 r— 1 r— 1 r— 1 

n n 

r— 1 r— 1 

Therefore if /(n) is divergent and i£(n) is an error term such that 

n n 

I>/(r) = yE/(r) + £7(n) 

r— 1 r— 1 

then linin-j.oo *£p^J(r) = ^' Therefore £7(n) = oQ^™ =1 f{r)). This completes 
the proof of the lemma. □ 

Lemma 2.5. If a n is uniformly distributed modulo one and g(n) is a monotonic 
and continuous divergent series of positive terms then, as n — >■ oo, 

n pi n 

E/KMr) = (1 + 0(1)) / /(x)dxE.9M- (8) 

Proof. The sum X)™=i f( a r)9( r ) wm depend on the sequence o r , n, the func- 
tion / and g, and may be also on some other hidden function say h which is 
independent of / and g. Hence without loss of generality we can express the 
above sum as 

n 

E /KMO = F(n)G(n)H(n) + E(n) (9) 

r=l 

where F(n) depends only on /, G(n) depends only on g and H(n) includes all 
terms which are not contributed by / and g; this may include constants and 
other terms (if any); and E{n) is an unknown error term. If suppose any one 
or more terms of the factorization F(n)G(n)H(n) did not exist then the corre- 
sponding value of F(n) or G(n) or H (n) will be unity. Later we shall see that 
H(n) = 1; but at this stage we do not know this. 

Our approach will be in two steps. In the first step, we choose the function 
g and impose no restriction on /. Since there is no restriction on / hence in the 
resulting sum, F(n) and H(n) will be the general form of F and H respectively 
while G(n) will be special case of G for the chosen function g. 

In the second step, we choose the function / and impose no restriction on g. 
In this case, G(n) and H(n) will be the general form of G and H respectively 
while F(n) will be a special case of F for the chosen function /. Since by defini- 
tion F, G and H are independent of each other, therefore after these two steps 
we will get the general forms of F(n) and G(n). 
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Step 1. Taking g(x) = 1 in [9] we get 

n 

J2f(ar)=F{n)G 1 (n)H(n)+E 1 {n) 



r=l 



where G\{n) and Ei(n) denote the function G(n) and E(n) respectively in the 
special case when g(x) — 1. From[T] Ei(n) — o(n). Hence, 

F(n)H(n) = ^ ( £ f(a r ) o(n)j . (10) 
Step 2. Taking f(x) — x in[9j we get 

n 

Y, a r9{r) = F 1 (n)G(n)H(n)+E 2 (n) 

r=l 

where Fi(n) and E2(n) denote the function F(n) and E(n) respectively in the 
special case when f(x) — x. From [7] E 2 (n) — oQ^" =1 g{r)). Also since a n is 
uniformly distributed modulo one, (1 — a n ) is also uniformly distributed modulo 
one and therefore, as n — ¥ oo, 



5^(1 - a r )g{r) = Fi(n)G(n)H(n) + o( £ g( 

r—1 r—1 

n n 

£ ff (r) = 2F 1 (n)G(n)H(n) + o(^. 9 (r) 



G(n)H(n) = \ > j g(r) - o[ > ^ ff (r) ) ] . (11) 



From 10 with 5(2;) = 1 and f(x) any arbitrary function; and 11 with 
f(x) = x, g(x) any arbitrary function; we have 

F {n)G{n)H{nf = ^ /K) = gM&L^) Z "SIU &)) } 

2Fi(n)Gi(n) 

Multiplying the terms in RHS of |12| and simplifying using Lemma 2.1, we 
obtain 



„i n 

F(n)G(n)H(n) 2 = (1 + o(l)) / /(a;)dx V) S (r). 



(13) 



Notice that dominant term of in the RHS of [13] equation can be completely 
evaluated if we know / and g. Therefore there is no hidden term H (n) or we 
can simply say that H(n) — 1. This completes the proof of the lemma. □ 
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3 The Master Theorem 



With the above lemmas, we are now ready to unify divergent series, summa- 
tion of arithmetic functions and uniform distribution modulo one into a single 
relation which we call the Master Theorem. 

Theorem 3.1. If a n is uniformly distributed modulo one and d n is positive and 
the following conditions are satisfied: 

(i) g is Riemann integrable in [Si,S n ] where S n = Yl r =i ^ r > 

(ii) d x g{S x ) is monotonic and continuous in [l,n], 

then for every function f Riemann integrable in [0, 1], 

y2d r g(S r )f(a r ) ~ / g(x)dx f(x)dx. (14) 
r=1 JSi Jo 

Proof. The proof follows directly from [8] by replacing g(r) with d r g(S r ) and 
then using [4] The replacement is justified because d r g(S r ) is a function of r 
and is independent of f(a r ). □ 

Notice all the lemmas of Section 2 are now unified under the Master Theo- 
rem. We shall now apply the Master Theorem to derive some of the existing as 
well as new results in different areas of mathematics. In the rest of this paper 
the symbols and notations introduced in the Master Theorem shall hold good. 
Also, we shall assume that the conditions mentioned in the Master Theorem are 
satisfied. 



Example 3.2. 



^oo^^^^l' (15) 
« r=l 



1 \ ■> e 



H, 



lim - > — 7-- — — — , (16) 

n^oon^r 1 + a 4 v ' 

r—l v ' ' 

l ™\Yn V = (ln2) 2 . (17) 

n->-oo n 2 ^— ' (Inn — lnr)(l + a r ) v ' y ' 



4 The Equivalence Theorem 

Although the Master Theorem unifies divergent series with uniform distribution 
modulo one, we have a more explicit unification between the two. We shall prove 
that the ratios of the partial sums of a divergent series replicate the behaviors 
of sequence uniformly distributed modulo one. With the same notations and 



definitions as in [141 we have 



Theorem 4.1. If d n is positive and divergent and limn._j.oo S Z 1 = 1 the 



r—l 
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Proof. Proceeding exactly as in the proof of |4j we can show that 

t^Kk)=Ly(i> y+0 ( Km+ Mk)\y (19) 

Substituting x = y/S n and making n — > oo, 19 reduces to 

it d rf(*cr) = S n [ f(x)dx + 0(1^/(1)1 + 1^/(0)1). (20) 
r _j K ^n / Jo 

Since S n is divergent, dividing both sides of[20]by S n , we obtain 

ton J_^ dr/ (|l) = / /(x)^. □ 

Thus we have explicitly unified divergent series with uniform distribution 
modulo one. It must be noted that above theorem cannot have a general ex- 
tension analogous to the Master Theorem because the sequence of the ratios of 
the partial sums are not necessarily uniformly distributed modulo one. 

Example 4.2. 

lim Y dr = (21) 



n— too 



5 Asymptotic equidistribution modulo one 

A standard application of definite integrals that is taught in elementary calculus 
in evaluating limit of sums using the rectangle method 



1 n rl 

n-voo n ^— ' \n/ J 



f(x)dx. (22) 



This formula is similar to[TJ If we take d r = 1 in Theorem 4.1 we get both 
[1] and [22] Thus the Equivale nee Theorem not only explains why these two 
formulas are similar but also shows that they are actually two different forms of 
the same principle. Also since [22] holds for all functions / Ricmann integrable 
in [0,1] it implies that the sequence 

1 2 n 

; 5 '"J 

n n n 

approaches uniform distribution modulo one as n — > oo. Inspired form this 
example, we introduce the concept of asymptotic equidistribution modulo one. 

Definition 5.1. A sequence of numbers b n is said to be asymptotically equidis- 
tributed modulo one if the sequence of ratios 

bn b n bn 

approach uniform distribution modulo one as n — > oo. 
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A consequence of this definition, the following corollary follows trivially from 
the Equivalence theorem 18 that the following corollary holds: 

Corollary 5.2. If b n is asymptotically equidistributed modulo OTIC (Mid S n — 
J2r=l br then, 

n^oo S n ^ V b„ / n-yoo n ' \b n J J Q 
r—1 r—1 u 

This naturally brings the question: Which are the sequences that are asymp- 
totically equidistributed modulo one. In this regard, we have the following 
lemma. 

Lemma 5.3. Ifb n is a sequence such that lini„_ i . 00 -jpf- = t for every t, < t < 
1, then b n is asymptotically equidistributed modulo one. 

Proof. Let < a < b < 1. If the condition mentioned in the statement of the 
lemma is true then, 

lim b[M] ~ b[at] =b-a. 
n— >oo b n 

Therefore n — > oo the probability of finding an integer r such that a < b r /b n < b 
approaches b— a and hence b n is asymptotically equidistributed modulo one. □ 

Corollary 5.4. If two sequences are asymptotically equidistributed modulo one 
then their linear combinations are also asymptotically equidistributed modulo 
one. 

Proof. Trivial. □ 

Using this criteria, we can easily prove that the sequence b n = n c is asymp- 
totically equidistributed modulo one only if c = 1. Similarly we can prove that 
the sequence of harmonic numbers H n is not asymptotically equidistributed 
modulo one. A very useful application of Lemma |5.3| is in the sequence of 
primes. We have: 

Theorem 5.5. The sequence of primes is asymptotically equidistributed modulo 
one. 

Proof. From the prime number theorem, limn-^ ^ = 1. Hence, if < t < 1 
then, 

lim P -M=t+ lim t * at =t. 

n->oo p n n->oo Inn 



Therefore by Lemma 5.3 the sequence of primes is asymptotically equidis- 
tributed modulo one. □ 

Similarly we show from the asymptotic expansion of the n th composite num- 
ber (See [4]) that the sequence of composite numbers, c„, is asymptotically 
equidistributed modulo one and consequently, we have: 
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Lemma 5.6. If p n and c n denote the sequence of prime numbers and the se- 
quence of composite numbers respectively and a, f3 and 7 are constants, such 
that ap n + /3c n + 771 ^ then, 

ii m -E/( ^ + f r+7r ) = fmdx. 

Proof. Since the sequence p n , c n and n all satisfy the condition of Lemma [5. 3| 
they are all asymptotically equidistributed modulo one and hence their linear 
combination will also be asymptotically equidistributed modulo one. □ 

Example 5.7. 

1 (P\ a ~ 1 ,„b-l ( cx \ c ' 



lim "TtE (-) ln "" - = -T^alnc) (23) 



n b+l ab+a 



^™ ^ r ab+a+6 

nH^ob J2r=i r a (l a + 2 a + . . . + r a ) b ~ ( a + 1/1 ( 25 ) 



1 n 2 i rl 2 i™2 

lim -E^h= iim -E-» £? h» = iim -£-^ = 5 ( 26 ) 



where T(u, v) denotes the incomplete gamma function in 23 and a > 1 in 24 



6 Summation over prime numbers 

By defining d r as a characteristic function that vanishes whenever r does not 
have a predefined characteristic, we can adapt the Master theorem to evaluate 
summation over prime numbers by defining d r — when r is composite. With 
this adaptation we have: 

Lemma 6.1. If p n is the n th prime number then, 

« rS Pr rl 

2]d Pr g(S Pr )f(a r ) ~/ g(x)dx f(x)dx. (27) 
r=1 Js 2 Jo 

Proof. The proof follows directly by replacing d r with d Pr in the Master Theo- 
rem [TU □ 

Lemma 6.2. //p„ is the n th prime number then, 

J2d Pr g(S Pr ) = / fl(tf)dy + 0(5(n)). (28) 

where 0(8(n) is defined analogous to the definition given m[4j 

Proof. The proof follows directly from [4| □ 
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Example 6.3. If 6{x) denotes the Chebyshev function of the first kind then 

^m =lax+O0) - (29) 

p<x yF > 

Proof. Form the Prime Number Theorem we have lim^^oo = 1. Hence, 
taking d Pr — lnp, and g(x) — ~ in [4J we have 

5^^ = ^) + 0(1)= lux + 0(1). 

p<x 



By actual computation, we have found that 

lim ( -^+\- — In x) 
x^oo V ^ 6(n) J 



converges to a constant value 8 sa 0.50904 . . .. The details of the calculation of 
the constant 9 are given in section 10. □ 

Example 6.4. Similarly we can show that ifip(x) denotes the Chebyshev func- 
tion of the second kind then 

^|g= In, + 0(1). (30) 

p v <x 

7 Consequences of Merten's Theorem 

One of the most celebrated result in the theory of prime numbers is the Merten's 
Theorem on the sum of the reciprocal of primes; 

V- = hxhix + M + o(z — Y (31) 
+^ p V In x ) 

p<x 

where M « 0.2614972... is the Meissel-Mertens constant(See [5], Theorem 
429) . In this section, we shall apply the Master Theorem on Merten's Theorem 
to derive several beautiful identities on prime numbers. Taking d Pr = l/p r in 
[28l we have 



Lemma 7.1. If f is monotonically decreasing then 



E -/(- + - + ... + -) ~ / " f(x)dx (32) 



<x Pr ^Pl Pi Pr> Jo.5 



In In x+M 



Proof. Taking cL = l/p r in 28 and simplifying using 31 we obtain the required 



result. □ 
Corollary 7.2. If c > 1 and m = [\np n ] then 

"\ r M+H r « r f+Qr r 1+M pi 

E — - — - E — /K) - -irr- / ( 33 ) 

where H r = 1 + I + . . . + 1 and O r = ^- + ^- + . . . + ^- . 
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Example 7.3. Taking c = e, f(x) — 1 and m — [\np n ] in 33 we obtain 
m e M+i+i+...+i ™ e -'+^ + ^+--+ih / 1 



E 



E 



Pr 



i? + o 



\np 7 



(34) 



r— 1 r— 1 

where H is a constant. Calculation of the constant H is given in Section 9. 



Example 7.4. Taking d Pr = 1/r + l/p r in 27 we obtain the identity; given in 
the introduction in section 1; that connects all the four fundamental constants 
e, 7r, 7 and M . 



hm — > - + — 

n->oo p n \ r p r / 1 + a r 



Tie 



7+M 



(35) 



where G r = 1 



I + ± + ± 

r pi p 2 



Example 7.5. Taking c = e and f{x) = 1 in 33 and comparing the result with 
Merten's theorem Y\ v<x ^1 — ~ fn~F> we obta-in 



g ' PI P2 



0(1). 



(36) 



The above identity is interesting because it connects a product involving the 
first n primes to a sum involving the first n primes. By actual calculation, we 
found that O-constant which we denote by p ~ 0.76774 .... 



8 A beautiful formula 

In this section, we shall demonstrate how our idea of unifying different area of 
mathematics can be used to yields beautiful results. As an example, we shall 
derive the intriguing formula that we have seen in the introduction in section 1. 



-4— E — f " + — ) ( tan 1 —) 2 ex P ( H r + — + — + ... + — 
m n 7r r Pr 7 ™ Pl P2 Pr 



3 7+M / G 7C(3)\ „/lnlnn 



4 

This formula brings together the elements from nine different topics of num- 
ber theory into a single beautiful result. We have the Apery's constant ({3), 
Catalan constant G, Euler-Mascheroni constant 7, exponential constant e, har- 
monic numner H n , Meissel-Merten constant M, n th prime number p n , imaginary 
part of the n th non-trivial zero of the Riemann zeta function j n and the funda- 
mental constants of a circle tt. 

Since 7„ ~ 7 n is asymptotically equidistributed modulo one. Taking 

11 7r (tan" 1 a;) 2 
d r = — I ,a r = — ,f(x) = 

r Pr In X 

in the Master Theorem, the result follows. 
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9 Generalized zeta function and Euler constants 



Throughout this paper, we have used summation of the type ^ , d r f(S r ) 
where S r — J2k=i dk ■ This suggests that many mathematical properties could 
be unearthed by studying series of this type. Although there are several gen- 
eralizations of the Riemann zeta series and the Euler constants, in this section, 
we present a new generalization of each of these two terminologies. 

Definition 9.1. The weighted zeta function is defined as 

oo , 

aw=Ev ( 38 ) 

n=l 

Thus Riemann zeta series (See [6]) is the simplest of the family of weighted 
zeta series, corresponding to the special case with unit weights d n = 1. The 
Riemann zeta series is related to prime numbers through the Euler product; 
it would be interesting to investigate the analogues of prime numbers for the 
weighted zeta series, if there are any. 

Definition 9.2. // d r f{S r ) is positive, strictly decreasing and divergent then 
the generalized Euler constants are defined as 

ld x J(x), g (x)=J^(22d r f(a r )g(S r )- f{x)dx J g(x)dx). (39) 

The classical Euler-Mascheroni, 7 w 0.577215665..., corresponds to the 
special case when d x = f(x) = 1 and g(x) = l/x. 

Example 9.3. 

ll/x 1 x a — l{ a ) = 1™ ( / — ~ ) ■ 

' ' n->oo V ^ — ' r a + 1 / 

With this definition, we have 7(0) = 7, 7(1) w 0.8225 .... 



10 Computations 

We have performed some computer calculation to determine the 0(1) terms in 
some of the formulas as well as to check accuracy of some of them. All programs 
were written in Intel© Fortran and run on the 64 bits AMD® Opteron 2700 
MHz processors. 



10.1 Constant of 1291 

We will present the data for determination of 0(1) term in: 

The results are presented in the Table II; the values of x are in powers of 10 
and powers of 2, the largest is x — 2 46 = 7.037. . . x 10 13 . We infer from this 
Table, that in[29]6> w 0.50904.... 
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TABLE I 



X 


In p 

2-^p<x 9(p) 


T, P<X w$ - Hx) 


10000000 


15.6266542966473 


0.5085586456 


16777216 


16.1442710909375 


0.5087387574 








10000000000000 


29.4426443525143 


0.5090381435 


17592186044416 


30.0075139149438 


0.5090379703 


35184372088832 


30.7006623741728 


0.5090392489 


70368744177664 


31.3938126082563 


0.5090423025 



Again values of x formed the progressions x = 10™ and x = 2™ and last value is 
x = 2 46 = 7.04 x 10 13 and it took 1 month CPU time to reach this value. 



10.2 Constant of 

It is not easy to estimate the big-0 constant appearing in[30]because the Cheby- 
shev function ip{x) for each x has to be calculated separately, i.e. it is not pos- 
sible to use the value ip{ x ) to calculate ip(y) for y > x. It is clearly seen from 
the following formula: 

i>ix) = E lo s(p) = E a h = E L lQ g(^)/ lo e(p)J lo s(p) (4°) 

p k <x n<x p<x 



where A(n) is of course the Mangoldt function. In the last form in (40) the 
summand evidently depends on x. Due to this obstacle we have calculated 

only for arithmetic progression x — nl000000, instead of geometric progressions 
in the previous cases. The obtained data is presented in the Table IV. 



TABLE II 



X 


ln (p) 

Z—tp<x ip(p) 


Yj V <x ipll) m ( x ) 


1000000 


13.5176919008455 


-0.2978186571188 


2000000 


14.2109984371948 


-0.2976593013294 


3000000 


14.6162538296483 


-0.2978690169841 


4000000 


14.9037288388882 


-0.2980760801960 








120000000 


18.3046100971201 


-0.2983922036263 


121000000 


18.3128899746922 


-0.2984111288688 


122000000 


18.3211283777276 


-0.2984032249699 


123000000 


18.3292829091801 


-0.2984120041566 


124000000 


18.3374008408318 


-0.2983912827375 


125000000 


18.3454350973536 


-0.2983891979130 


126000000 


18.3534102043258 


-0.2983822605899 
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Above results suggest that in [30] 



E 

p v <x 



lnp 

ip(p) 



Ins + 0(1) 



there will be no equality, as LHS contains more terms and will be larger than: 

lnp 



E 

p<x 



1p(p) 



Indeed the computer has produced data: 



TABLE III 



X 


v - ^ In p 

2-^p v <x ip(p) 




1000000 


50.9115451893350 


37.0960346313707 


2000000 


54.6640936160706 


40.1554358775464 


3000000 


56.6781935153063 


41.7640706686740 


4000000 


57.2884765277029 


42.0866716086188 


5000000 


59.5704206070202 


44.1454721366218 


6000000 


60.2669838560583 


44.6597138288660 


7000000 


60.6473757815097 


44.8859550744901 


8000000 


61.0151186141244 


45.1201665144803 


9000000 


62.1874985975655 


46.1747634622650 


10000000 


62.7528494571091 


46.6347538061507 








26000000 


66.9185401883590 


49.8449330923732 


27000000 


66.9728912056416 


49.8615437816730 


28000000 


67.0256045468440 


49.8778894787045 


29000000 


67.2571282847139 


50.0743218967631 


30000000 


67.3132959212713 


50.0965879816448 


31000000 


67.3999358663197 


50.1504381038703 



10.3 Constant of 



We tried to determine the constant hidden in O(l) in 34 The results are not 
conclusive as the difference 

M+i+|+...+i ™ e "V+ST+ii+~+^ 



r i Pr 

r—1 r—1 

displays fluctuations between 5 and 7, see Table III 

TABLE IV 
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X 


ro e Af + l+i+... + l 

Z-^r—l r 


T+-L + J- + ...+ J- 
e pi p 2 Pr 


ClllieieilCe 


10000000 


41.0329563695293 


oa oo c c n /i o a >70 c no 
34.385594^4/8502 


6. 647362 Iz 16792 


16/ i (lib 


41.0329563695293 


or r o nr - fVl r 7 n i c on 

35. 5826945591639 


5.4502618103654 


OOCC/1/IOO 

33554432 


AO AA A 7AO A OO O C OT 

43.414 /08433352 / 


o^7 1 o£; or\o i i ooo i n 
3/ .1863081123819 


a ooo lAA'jnnA7AO 

6. 2284003209 /08 


6/ 108864 


45.7926459440599 


OO 7/On C ^700£27/00 70 

38. (895 / 286 1 2,6(6 


^7 AAOH70A7(;0007 

/.0030/30/682z7 










100000000000 


62.3588796562188 


55.6923315736851 


6.6665480825336 


137438953472 


62.3588796562188 


56.4280087636232 


5.9308708925955 


274877906944 


64.7168951615670 


58.0315248369373 


6.6853703246296 


549755813888 


67.0732528789821 


59.6350423017193 


7.4382105772628 


1000000000000 


67.0732528789821 


61.0190970811536 


6.0541557978285 


1099511627776 


67.0732528789821 


61.2385589160106 


5.8346939629716 


2199023255552 


69.4280715529841 


62.8420759971757 


6.5859955558084 


4398046511104 


71.7814576116916 


64.4455941070075 


7.3358635046840 


8796093022208 


71.7814576116916 


66.0491118126049 


5.7323457990866 


10000000000000 


71.7814576116916 


66.3458668229100 


5.4355907887816 


17592186044416 


74.1335068131232 


67.6526285355308 


6.4808782775924 


35184372088832 


76.4843056253357 


69.2561434465050 


7.2281621788307 



10.4 Verification of 1351 

We shall discuss the formula 

lim — V ( - + — ) — - = = 1.8169017889 . . . , 

n-^-oo p \r V r J 1 + dr 4 

r—1 

, 1 11 1 

G r = 1 + - + ... H 1 h • • • H 

2 r p 2 Pr 

which connects the fundamental mathematical constants tt, e, 7, ill and contains 
also sequence a r uniformly distributed on the interval (0,1). The numbers a r 
were generated from the sequence 

Z r +i = AZ r mod C, a r = — 



where A = 1203248318 and C = 2 31 - 1 (Mersenne prime). It is known, 
that period of Z n is C and the sequence is uniformly distributed over (0, 1) 
[9j. The period of a r was closed during calculations many times as there are 
346065536839 > 2 38 primes up to 10 13 , but we believe it is not an obstacle in 
this problem. The results are presented in the Table I; the values of x form the 
progression 10™ and 2™ — the last value is for x — 2 46 = 7.037 . . . x 10 13 . 



TABLE V 
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11 


y n (i + i ) e G \ 

Z-^r—l \ r p r J l-\-a r 2 


Z_/r=l y r T p ry / 1+a,. 2 4 


1.000000000 X 1U 


1 m a a*7 a o i r\ 
1.9144/431U 


n nn7C/i ^nco 
0.09/544058 


4.z949b/29i xiu 


1.9U f4zlz78 


0.09049102b 


i nnAnnnnnnn . . mlU 

l.UOUUOOUUOO xll) 


1.9U3( /3984 


n no^O/io^oi 
0.086843 /3z 


1 TITf^Q^'fAl A O . . 1 nil] 

i. ( i /y<5byi4o xiu 


i.yoi5/882y 


n no^cj/ioc'7'7 
0.084D48577 








1.000000000000 xlO 12 


1.888054174 


0.071123922 


1.0yy51162768y xlO 12 


1.887789938 


0.070859685 


2.199023255531 xlO 12 


1.885918614 


0.068988361 


4.3y80465110y3 xlO 12 


1.884144386 


0.067214134 


8.7y6093022151 xlO 12 


1.882458766 


0.065528514 


1.0000000000000 xlO 13 


1.882156156 


0.065225904 


1.7592186044399 xlO 13 


1.880856104 


0.063925852 


3.5184372088777 xlO 13 


1.879329708 


0.062399456 


7.0368744177643 xlO 13 


1.877873176 


0.060942923 



10.5 Constant off 



Writing 36 in the form L(n) = e'- M+0 ^R{n) + 0(1), where: 



and R(n) is RHS: 



(41) 



(42) 



\k=l 



we see, that p can be found by plotting points (L(n),R(n)) for some values of 
n and fitting the straight line to these points by the least square method. We 
computed the values of (L(n), R(n)) on the computer and stored to the file at 
n = 10 6 , 2 24 , 2 25 , . . . , 10 13 , 2 46 = 7.037... x 10 13 . The figure 1 presents obtained 
points. 



Fitting the straight line to the last 4 points (in fact two points determine line!) 
gives L = 1.371247673i? + 0.7676022888, while e 7 ~ M = 1.3712441303, thus 6 
digits of the correct value are reproduced. Assuming that for p also 6 digits are 
correct we have p = 0.76770. Inspection of the Sloane's Online Encyclopedia of 
Integer Sequences does not indicate what this number can be in terms of other 
mathematical constants. 
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10 20 30 40 

R(n) 

11 Future works 

The scope for future works should focus on extending the master theorem to 
other areas including the complex domain. In particular since uniform distribu- 
tion is a special case of low discrepancy sequences, we would want to extend the 
results of this paper to low discrepancy sequences in higher dimensions. Specific 
area of interest would be the Quasi-Monte Carlo method and particularly the 
Koksma-Hlawka Inequality (See [7]). 
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